Tandem Pairing in Heavy Fermion Superconductors 



Rebecca Flint and Piers Coleman 

Center for Materials Theory, Rutgers University, Piscataway, NJ 08855, U.S.A. 

We examine the internal structure of the heavy fermion condensate, showing that it necessarily 
involves a d-wave pair of quasiparticles on neighboring lattice sites, condensed in tandem with a 
composite pair of electrons bound to a local moment, within a single unit cell. These two com- 
ponents draw upon the antiferromagnetic and Kondo interactions to cooperatively enhance the 
superconducting transition temperature. The tandem condensate is electrostatically active, with a 
small electric quadrupole moment coupling to strain that is predicted to lead to a superconducting 
shift in the NQR frequency. 
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In many strongly interacting materials, quasiparticles 
are ill-formed at the superconducting transition, giving 
the Cooper pair a non-trivial internal structure. The 115 
family of heavy fermion superconductors [T]-[4] provide an 
extreme example of this phenomenon, where quasiparti- 
cle formation, through the screening of local moments by 
electrons, coincides with the onset of superconductivity. 

The 115 family has long attracted great interest for 
the remarkable rise of the superconducting transition 
temperature from T C =0.2K in Celn3 under pressure [5] 
to 2.3K in CeCoIn 5 [l-3 and then up to 18. 5K in 
PuCoGa5[4 . While the abundance of magnetism in the 
phase diagram has led to a consensus that spin fluctu- 
ations drive the superconductivity in the cerium com- 
pounds [5HS] , the presence of unquenched local moments 
at T c is difficult to explain within this picture. In a typ- 
ical spin-fluctuation mediated heavy fermion supercon- 
ductor, the local moments quench to form a Pauli param- 
agnet (x(T) ~ xo) we U before the development of super- 
conductivity. Yet NpPd 5 Al 2 [10] and Ce{Co,Ir}In 5 [2j [11] 
exhibit a Curie- Weiss susceptibility, x(T) ~ 1/(T+Tcw) 
down to T c . Moreover, the highest transition tempera- 
tures are found in the actinide 115s, which show no signs 
of magnetism. 

These observations led us to recently propose [12] that 
the actinide 115s are composite pair superconductors [T3] . 
where the heavy Cooper pair forms by combining two 
electrons in two orthogonal Kondo channels with a spin 
flip to form a composite pair, Ac = (N\c\^C2^S+\N + 

2), where c[ 2 create electrons in two orthogonal Kondo 
screening channels p~2j[14]. However, composite pairing 
alone cannot account for the importance of magnetism 
in the Ce 115 phase diagram. 

We are led by these conflicting observations to propose 
a model for the 115 materials where the composite and 
magnetic mechanisms work in tandem to drive supercon- 
ductivity. Composite pairing originates from two channel 
Kondo impurities, while magnetic pairing emerges from 
antiferromagnetically coupled Kondo impurities. These 
two systems are equivalent at criticality in the dilute 
limit [15,, and we argue that this connection persists to 



the lattice superconducting state that conceals a common 
quantum critical point (QCP)[T6], 
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FIG. 1: (Color online) A tandem pair contains a superpo- 
sition of magnetic and composite pairing, both with d-wave 
symmetry. The magnetic pair (left) contains neighboring f- 
electrons, while the composite pair (right) combines a spin 
flip and two conduction electrons. The unit cell is denoted by 
dotted lines, with dots indicating the local moment sites. 

To expose the interplay between magnetic and com- 
posite pairing, we examine the internal structure of a 
heavy fermion pair. In a Kondo lattice, the heavy quasi- 
particles are a linear combination aj^ = u^c^ + v^f^, 
where and p create conduction and localized elec- 
trons, respectively [1 7 . The wavefunction is 



|*>=P G exp(AT)|0>, 



(1) 



where A^ = ^ k A^a^al^) creates a d-wave pair of 
quasiparticles and Pq is the Gutz wilier projection oper- 
ator restricting the number of f-electrons to one. Act- 
ing the Gutzwiller projector on the f-electron reveals 
its internal structure as a composite between a conduc- 
tion electron and a spin flip at a given site j, PgIj^ ~ 
(c^S+)Pg- The pairing field contains three terms 
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The diagonal terms, with A£ = u^A^ and A^ = ^ A k 
create f- and conduction electron pairs. A d-wave pair of 
f-electrons is an inter-site operator, taking the form 



*L = ^A M (R ii ) (4^S i .)(c] i S j+ ) 



(3) 



outside the Gutzwiller projection. However, if we expand 
the off-diagonal terms in real space, 



*t, = £A c (Ry) 



(4) 



where A C (R) = Xlk(' u k^kAk)e 2k ' R , we find a composite 
pair formed between a triplet pair of conduction electrons 
and a spin flip [T2HT4] . Unlike its diagonal counterparts, 
which are necessarily inter-site, composite pairs are com- 
pact objects formed from pairs of orthogonal Wannier 
states surrounding a single local moment (Fig. [T]). 

Magnetic interactions favor the inter-site component 
of the pairing, while the two-channel Kondo effect fa- 
vors the composite intra-site component. However, both 
components will always be present in the superconduct- 
ing Kondo lattice. If the product of the Kondo screening 
channels has a d-wave symmetry, the composite and mag- 
netic order parameters necessarily couple linearly to one 
another, a process that enhances the transition tempera- 
ture over a large region of the phase diagram, providing 
a natural explanation for both the actinide and Ce 115s. 

To treat these two pairing mechanisms simultaneously, 
we introduce the two channel Kondo-Heisenberg model, 

H = H C + H K1 + H K2 + H M (5) 

and solve it in the symplectic-A^ limit [12], where 
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where Sj is the local moment on site j, and ipjr is the 
Wannier state representing a conduction electron on site 
j with symmetry T, 



^Ta = X^ Fka6Ck6eikRi ' 



(8) 



where the form factor &rkab is only diagonal in the spin 
indices in the absence of spin-orbit. Microscopically, the 
two orthogonal Kondo channels, Jr arise from virtual 
fluctuations from the ground state doublet to excited sin- 
glets, where the two channels correspond to adding and 
removing an electron, respectively. The Ce 4/ 1 state is 
split by tetragonal symmetry into three Kramer's dou- 
blets, where is the ground state doublet [T8l [19], so 
we may summarize the virtual valence fluctuations with: 
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Requiring the composite pairing to resonate with the d- 
wave magnetic pairing [20] uniquely selects (8)T 6 as the 
lowest doubly occupied state, as this combination leads 
to d-wave composite pairing[12]. A simplified two dimen- 
sional model is sufficient to illustrate the basic physics, 
where the d-wave composite pair now comes from the 
combination of s-wave hybridization in channel one and 
d-wave hybridization in channel two [21] [22]. The mag- 
netism is included as an explicit RKKY interaction, Jh 
between neighboring local moments (ij), generated by 
integrating out electron in bands far from the Fermi 
surface. Treating the magnetism as a Heisenberg term 
leads to a two band version of resonating valence bond 
(RVB) superconductivity [23], where the local moments 
form valence bonds which "escape" into the conduction 
sea through the Kondo hybridization to form charged, 
mobile Cooper pairs [24]. 

To solve this model, we use a fermionic spin represen- 
tation, Sj = fj&Nfj'i symplectic-A^ maintains the time- 
reversal properties of SU (2) in the large N limit by using 
the symplectic generalization of the Pauli matrices gn to 
construct the spin Hamiltonians [T2] . 
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where e is the large N generalization of ia2 • Each quartic 
term can be decoupled by a Hubbard-Stratonovich field, 
leading to normal, Vp and anomalous, Ar hybridization 
in each Kondo channel and particle-hole, hij and pairing, 

The SU(2) gauge sym- 
->• uf + ve} ft is used to 
eliminate Ai. The lowest energy solutions contain only 
pairing fields in the magnetic and second Kondo chan- 
nels, giving rise to three Hubbard-Stratonovich fields, 
Vi, A2 and A#, where Ah is d-wave in the plane, so 
An(cosk x — cos k y ). Using the Nambu no- 
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the mean field Hamiltonian can be 
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where A is the Lagrange multiplier enforcing the con- 
straint nf = 1. The mean field Hamiltonian can be diag- 
onalized analytically. Upon minimizing the free energy, 
we obtain four equations for A,Vi,A 2 , and A#. Solv- 
ing these numerically, and searching the full parameter 
space of J2/J1, Jh/ Ji and T to find both first and sec- 
ond order phase transitions, we find four distinct phases: 
a light Fermi liquid with free local moments when all pa- 
rameters are zero, at high temperatures; a heavy Fermi 
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liquid when either V\ or A2 are finite, with symmetry 
T, below Tkt] a spin liquid state decoupled from a light 
Fermi liquid when Ah is finite, below Tsl] and a tandem 
superconducting ground state with Vi, A2 and Ah all fi- 
nite, below T c , as shown in Fig. [2] There is no long range 
magnetic order due to our fermionic spin representation. 




FIG. 2: (Color online) The superconducting transition tem- 
perature as the amounts of magnetic, Jh and second channel, 
J 2 couplings are varied. A slice at T = Tki shows the regions 
of the spin liquid and Fermi liquids, and the orange ellipse 
is a path illustrating how materials could tune the relative 
coupling strengths (see Fig. [3]). The phase diagram was cal- 
culated in a simple two dimensional model with channel one 
s-wave and channel two d-wave (n c = .75). The transition is 
first order for Jh/ Ji > 4, but otherwise second order. 

Experimentally, CeMIn 5 can be continuously tuned 
from M = Co to Rh to Ir[3 . While CeRhIn 5 is a 
canonical example of a magnetically paired superconduc- 
tor, where moderate pressure reveals a superconduct- 
ing dome as the Neel temperature vanishes pQ, further 
pressure [25] or Ir doping on the Rh site [3! leads to a sec- 
ond dome, where spin fluctuations are weaker [26 . We 
assume that the changing chemical pressure varies the 
relative strengths of the Kondo and RKKY couplings, 
so that doping traces out a path through the phase dia- 
gram like the one in Fig. [3J chosen for its similarities to 
CeMIn 5 . While different paths may lead to one, two or 
three superconducting domes, by maintaining the same 
Fermi liquid symmetry throughout (Tki > Tk2), we are 
restricted to one (magnetic only) or two (magnetic and 
tandem) domes. 

A qualitative understanding of this tandem pairing can 
be obtained within a simple Landau expansion. For T ~ 
T c <C Tki-, $ = A 2 and \I/ = Ah will be small, and the 
free energy can be expressed as 
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FIG. 3: (Color online) A possible experimental path through 
the phase diagram in Fig [2] chosen for its similarity to the Ce 
115 doping phase diagram [3]. The transition temperatures 
for superconductivity (T c in solid blue), spin liquid (Tsl in 
dotted red), and Fermi liquid (Tki in dashed orange and Tki 
in dot-dashed white) are plotted for comparison. Tempera- 
tures are scaled by Tki, which may itself vary as one moves 
around the phase diagram [27]. While we always find a super- 
conducting ground state, due to our choice of a fermionic spin 
representation, real materials will have an antiferromagnetic 
ground state for Tsl /Tki sufficiently large. 



coupling of the two order parameters, 7 = dF/dA2dAn 
is always nonzero in the heavy Fermi liquid, leading to 
an enhancement of the transition temperature, 
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<^i,2, Pi, 2,1 an d 7 are all functions of A and V\ and can 
be calculated exactly in the mean field limit. The linear 



For /?i/?2 > fif, the two order parameters are only weakly 
repulsive, leading to smooth crossovers from magnetic to 
composite pairing under the superconducting dome [28 . 

While the development of conventional superconduc- 
tivity does not change the underlying charge distribu- 
tion, tandem pairing is electrostatically active, as com- 
posite pairing redistributes charge, leading to an electric 
quadrupole moment. The transition temperature of the 
115 superconductors is known to increase linearly with 
the lattice c/a ratio [29], conventionally attributed to de- 
creasing dimensionality. Our theory suggests an alter- 
native interpretation: in a condensate with a quadrupole 
moment, Q zz oc ^3/^, which couples linearly to the tetrag- 
onal strain, AF oc —Q zz Ut e t, the second term in the Lan- 
dau free energy (12) becomes a 2 [T — (T c2 + \u te t)]^ci 
naturally accounting for the linear increase in T c . The 
development of a condensate quadrupole moment should 
be also detectable as a shift of the nuclear quadrupole res- 
onance (NQR) frequency at the nuclei of the surrounding 
ions. 



4 



The link between f-electron valence and the Kondo 
effect is well established [30], but tandem pairing intro- 
duces a new element to this relationship. Changes in the 
charge distribution around the Kondo ion can be read off 
from its coupling to the changes in the chemical poten- 
tial, Ap(x) — \e\SH I 'Sja(x). The sensitivity of the Kondo 
couplings to /i is obtained from a Schrieffer- Wolff trans- 
formation of a two-channel Anderson model, which gives 
J' 1 = AE F /V* . Here, Vp. o are the bare hybridiza- 
tions and A£r are the charge excitation energies. With 
a shift in /x -> /x + 5p(x), SJ^ 1 = ±\^r(x)\ 2 Sfi(x)/V 2 . 
The sign is positive for J\ and negative for J 2 because 
they involve fluctuations to the empty and doubly occu- 
pied states, respectively: /° ^ f 1 ?=± f 2 . Differentiating 



(11) with respect to 5/j,(x), the change in p(x) will be: 



Ap(x) = \e\ 



A. 
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(14) 

For equal channel strengths, the total charge is constant, 
and the f-ion will develop equal hole densities in and 
electron densities in T^, leading to a positive change in 
the electric field gradient, dE z /dz oc (T c — T) > at the 
in-plane In site that will appear as a shift in the NQR 
frequencies growing abruptly below T c (see Figure P}. 
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FIG. 4: (Color online) Predicted NQR frequency shift, 
Avnqr in CeMIns. The inset shows the relative locations of 
the indiums in-, In(l) and out-of-plane, In(2). Avnqr mea- 
sures the change in the electric field gradient (EFG) due to 
the onset of superconductivity. For equal channel strengths, 
the total charge of the f-ion remains unity, but the increas- 
ing occupations of the empty and doubly occupied sites cause 
holes to build up with symmetry (orange) and electrons 
with symmetry T6 (blue). The change in charge distribution 
and resulting electric fields are shown above in a slice along 
the [110] direction (dashed line in the inset). The positive 
EFG, dE z /dz at the In(l) site will lead to a sharp positive 
shift in vnqr starting at T c . 

The f-electron valence should also contain a small 
superconducting shift, observable with core-level X-ray 



^ 2 C oc (T c - T), as ^ c oc A 2 when J x > J 2 . While the 
development of Kondo screening leads to a gradual va- 
lence decrease through Tk, as it is a crossover scale, the 
development of superconductivity is a phase transition, 
leading to a sharp mean- field increase. Observation of 
sharp shifts at T c in either the NQR frequency or the va- 
lence would constitute an unambiguous confirmation of 
the electrostatically active tandem condensate. 
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